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Fig. 1 Equilibrium profiles of the magnetic (dased lines) 
and velocity fields (solid line). 



The dynamics of a global reconnection in the presence of a poloidal shear flow which is located in 
between magnetic islands is investigated. Different linear regimes are identified according to the value of 
the resistivity and the distance between the low-order resonant surfaces. It is found that the presence of 
a small shear flow affects and significantly delays the global reconnection processes. It is shown that this 
delay is linked to a breaking of symmetry imposed by the existence of the shear flow and the generation 
of a mean poloidal flow in the resistive layers. 
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1. Introduction 

In tokamaks, the internal transport barrier is 
observed in reversed magnetic shear configurations 
where two low-order rational surfaces exist [Tj . In such 
configurations, toroidal and poloidal flows, as well as 
temperature and density gradients can coexist close 
to the plasma core [2]. Indeed, there is evidence that 
such zonal flows exist and that the shear flow is lo- 
calized in between resonant surfaces where the double 
tearing instability can grow. The observed maxima of 
the flow velocity are weak compared with the Alfven 
velocity va- Nevertheless, such flows can be impor- 
tant with respect to the turbulence level because the 
associated radial electric field can be linked to the for- 
mation of an ITB[3]. In fact, such configurations, al- 
lowing eventually the generation of a strong internal 
barrier, are part of the ITER scenarios for advanced 
confinement [I] . 

Previous works have investigated the influence of 
a Bickley jet on a tearing instability [51 [S] and on a 
double tearing instability [7]. More recently, Bierwage 
et al [H] have studied the influence of the core ro- 
tation amplitude - for a given sheared flow profile 
- on the stability of magnetohydrodynamic (MHD) 
and Kelvin-Hclmholtz instabilities (KHI) . These stud- 
ies are mainly linear and a cylindrical geometry is 
used. They have shown that KHI can develop at large 
poloidal mode numbers and are strongly enhanced 
when the core rotation passes a critical value of the or- 
der of 1Q~ 2 va in the case where the distance between 
the resonant surfaces is small. They also found that 
below this critical value the growth rates of both the 
MHD and KH mode decrease with the amplitude of 
the core velocity. 



In this paper, we focus on the case where the shear 
flow is in between the two tearing instabilities. Within 
the framework of MHD, a 2D slab geometry is used. 
It is found that the nature of the dynamics depends 
drastically on the distance between the low-order reso- 
nant surfaces. The impact of the flow on the nonlinear 
process leading to a global magnetic reconnection of 
the system is studied. The origin and the role of the 
dynamically generated m = poloidal flow turbulence 
are investigated as well as the dynamics at small scales 
once the reconnection has taken place. 

2. Model equations 

We use a two field model corresponding to a re- 
duced (MHD) description of the fluid equations that 
provides a minimal framework to study the impact of a 
radially sheared flow 4>o{x) on a double tearing mode. 
The model consists of a set of two coupled equations 
for the fluctuations of the electrostatic potential 4> and 
magnetic flux if). The equilibrium magnetic field in the 
z— direction is given by the constant Bq z and in the 



poloidal by Bo(x) = ip' (x)y. The equilibrium current 
is therefore jo = V'o ( x )- The time evolution of the two 
fields is described by 

8 t UJ+ [<p + (f>Q,UJ + UJ Q ] = [tp + l/jQ, j + j ] 

+vV\u , (1) 

d t iP+ [0 + 00,^ + ^0] = m , (2) 

where r/ is the resistivity, v the viscosity, u> = 
the vorticity, and j — S/ 2 j_ip t ne current density fluctu- 
ation. The equations are normalized by ta = L±/va 
for the time, L±Bq z for ?/>, and Lj_va for <f> where ta 
is the Alfven time and L± is a magnetic shear length. 

The poloidal equilibrium magnetic field Bq(x) is 
chosen as 

Bo = tanh f ?—IlL) _ tanh ( ^^L\ (3) 
V as ' \ as J 

where the parameter as = 0.5 controls the width of 
the profile, x G [— L x /2, L x /2] and {x\ ri xir — —x\ r ) 
are the locations of the resonant surfaces where double 
tearing instabilities develop. The poloidal equilibrium 
shear flow is given by 

wo = A v tanh(a;/a. u ) , (4) 

where A v = 0.03 and a v — 0.2. The resulting pro- 
files are shown in fig. fl}. The dashed lines rep- 
resent four typical equilibrium magnetic profiles we 
have used, corresponding to Sx — xi r — x\ r — 
{ir/2, 37r/4, 7r, 57r/4} (cyan, red, green and blue re- 
spectively). The purple curve indicates Vq. Equa- 
tions (1~2) are solved numerically using a finite differ- 
ence scheme in the x-direction, including an Arakawa 
algorithm for an accurate conservation of the Pois- 
son brackets [.,.] and a pseudo-spectral method in 
the y— direction, including an appropriate de-aliasing 
scheme[9l 110) . The number of grid points in the x 
direction range from N x — 256 to N x — 2048, with 
a typical spatial resolution of the order dx = 5.10~ 3 . 
L x and L y are the box size in the x and y directions 
respectively. We find that these different profiles give 
rise to very different dynamics. For instance, when 
Sx = 57r/4, no global reconnection occurs and complex 
nonlinear dynamics linked to the poloidal rotation of 
the islands are observed. When Sx = n/2, islands do 
not rotate and there is a global reconnection process 
in between the resonant surfaces. In this paper, we 
focus mainly on the latter situation. 

3. Coexistence of linear instabilities 

It is well known that a constant magnetic field 
stabilizes a shear flow instability when the initial flow 
is a vortex sheet [TT] , as far as the Alfven velocity ex- 
ceeds the amplitude of the vortex sheet. In [7], it is 
shown numerically that when the flow is forced with a 
velocity field v(x) = A(l/ cosh(a;/a) — l)y, KHI coex- 
ist with the magnetic tearing whenever a < 1. In the 
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Fig. 2 Scheme for a calculation: Without shear flow (blue 
line), 7(log(?7)) is a line and a the associated slope. 
With a shear flow, two different branches exist. For 
the low resistivity branch, KH is unstable (green 
line) and for the high resistivity one, KH is stable 
(red line). 
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Fig. 3 a as a function of the distance between rational 
surfaces Sx. 

latter work, the double tearing mode in the presence 
of a shear flow in cases where no KHI develops has 
also been studied, while in [5] this has been done in a 
case where KHIs are present. However, the conditions 
at which KHI and DTM coexist are not yet clear. 

It is instructive to focus on the role of both the 
distance between the resonances Sx and the resistivity 
for a fixed amplitude of the flow A v . In the absence of 
shear flow, the growth rate of the resistive DTM scales 
with the resistivity like r/ a with a = 3/5 in the limit of 
very ideally stable MHD conditions and a = 1/3 when 
the system enters an ideally marginally stable regime. 
Ideally Marginally stable means that the free mag- 
netic energy of the ideal mode A# is zeroJTJ [TJ]. In 
the limit of closed resonant surfaces, the latter can be 
linked to the tearing instability parameter A'. Xh de- 
pends on the magnetic equilibrium profile and there- 
fore on Sx. Fig. [3] shows the dependance of growth 
rate on the resistivity in the case with/ without shear 
flow. Changing Sx and rj |13j . we plot the power law 
index a{Sx) in fig. ([3]) . As expected, without shear 
flow, the power law parameter a ranges in between the 
two limit cases 1/3 and 3/5 (blue curve). When the 
distance between two rational surfaces is large enough, 
we converge to the standard tearing law. Conversely 
when Sx becomes of the order of a typical magnetic 
shear length L±, wc enter into a full DTM regime. 



A global magnetic reconnection occurs nonlinearly in 
such cases. 

When a weak shear flow is added, we have to dis- 
criminate between two regimes. A regime in which 
the KHI is unstable appears, depending on the nu- 
merical value of the resistivity. This instability is ra- 
dially localized in the vicinity of the layer where the 
velocity shear is maximal x ~ 0. Typically, it has 
a high poloidal mode number but we observe that 
it modifies the growth of the tearing m = 1 mode 
(fc m = to 2w/L y ). In fact, this not surprising because 
the DTM regime is linked to a strong radial coupling 
between the two magnetic surfaces and the presence 
of vortices in between them modifies the nature of the 
interaction, and a priori should weaken it. In fig. ([3]), 
the high and low resistivity cases, the KH stable (red 
curve) and KH unstable regimes (dark green curve) 
are shown. We observe that, for Sx < 2, the presence 
of a KH instability in between the resonances amplifies 
the power law index strongly such that it even exceeds 
the asymptotic regime a = 3/5. On the other hand, 
we observe of strong decrease of the power law index 
for 2 < Sx < 4 in the KH stable regime, wich roughly 
follows the case without shear flow, but is even lower. 

The reason is that the imposed shear flow shown 
in fig. (fl]) ideally leads to a poloidal rotation of the 
island in both directions. The global poloidal rota- 
tion of the plasma modifies the energy balance of the 
system and in fact reduces the growth rate of the mag- 
netic instability. The viscosity, the resistivity, and 
also the nature of the mode can prohibit such an ideal 
scenario. In fact, the plasma starts to rotate around 
Sx ~ 2, but when the resonances are closer together, 
given the weakness of the amplitude of flow A v and 
the strong interaction between the growing islands, 
the plasma poloidal rotation is locked. 



4. Structure of the m=l DTM and gen- 
eration of a mean poloidal flow 

In this section we focus on the nonlinear evolu- 
tion of the m=l DTM. Fig. 2] shows the time evolu- 
tion of the kinetic and magnetic energy of the m = 
{0,1,2,3,4,12} modes for the case with Sx = n/2, 
-q = 10~ 3 , /i = 10~ 4 , L x = 2-7T, and L y — 2ir, where 
KH is linearly stable and only the m=l DTM is un- 
stable. The m > 1 modes are stable in the linear 
phase (t < 160) and are destabilized in the quasilin- 
ear phase, where the exponential growth of the modes 
satisfies the relation -f m — mji (t < 400) . Contrary 
to the vq — case where the to = mode is not gener- 
ated in the linear and quasilinear phases, the relation 
7o = 271 is held in those phases. We will see that the 
generation of the to = mode modifies the dynam- 
ics of the system leading to the global reconnection in 
between the resonances. In order to understand the 
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Fig. 4 Energy of the poloidal modes versus time. 5x 
tt/2. 
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Fig. 5 Snapshots of (Left) and ip (Right) at t = 220ta- 
Sx = it/2. 



origin of this mode, it is instructive to observe the 
snapshots of the electrostatic potential and the mag- 
netic flux in the quasilinear phase. They are shown 
in fig. ([5]). Note that the structures are not sym- 
metric with respect to a reflection at x — and/or 
y = constant as discussed in the previous section. De- 
spite an imposed poloidal shear flow, the islands do 
not rotate poloidally. However, a symmetry is clearly 
broken when compared with the vq — case. Indeed, 
in the latter case, the m = 1 mode satisfies a symme- 
try (f>i(x, -y, t) = -<pi(x, +y, t) where 



(f>i(x,t) cos(a%(x) + k%y) (5) 



is the m = 1 mode, and <pi and a,\ are the amplitude 
and phase of the mode respectively. In other words, 
when vq — 0, a% is constant because of the symme- 
try. This is clearly false when vo ^ as is seen in 
fig. OH In fact, a straightforward calculation shows 
that whenever a% is constant, the projection of the 
Poisson brackets [0i,a->i] on to = mode is zero, for 
the same reason the projection of will also be 

zero when vq = 0. From the breaking of poloidal par- 
ity as soon as vq 0, it follows that the generation of 
the m = mode occurs at, of course, a growth rate 
271 through both the Reynolds and Maxwell stresses. 




Fig. 6 Total magnetic flux just before and at the end of 
the global reconnection process. Sx — tt/2. (Left) 
t = U0t a . (Right) t = 470t a . 
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Fig. 7 Sx — 7r/2. Perturbed electrostatic potential just 
before the global reconnection. (Left) No shear 
flow. (Right) With the shear flow 



5. Impact of the shear flow on the 
global reconnection process 

When the islands are sufficiently closed together, 
the global reconnection occurs regardless of a shear 
flow, as is seen in fig. [5] where the nonlinear interac- 
tion of the two magnetic islands deforms the DTM 
structure (t = 440ta) and forces them to vanish pro- 
gressively (f = 470ta)- ft is interesting to check if 
this process is linked to the generation of the mean 
velocity flow fluctuation vq = VQ at (x,t) — vq(x) in the 
quasilinear phase of the m = mode. A close exam- 
ination of fig. studied together with snapshots of 
the magnetic flux shows a correlation. Indeed, the tp\ 
component of the m = 1 mode starts to move radi- 
ally when the kinetic energy of the m = mode, wich 
is linked to vq, becomes more important than that of 
the m = 1 mode (step 1), in the Sx — tt/2 case at 
t ~ 395. The global reconnection process can then 
develop and occur in the time interval t 6 [440,480]. 
In fact, when the magnetic energy of the m — mode 
crosses that of the m = 1 mode at t ~ 448 (step 

2) , the remaining islands are no longer topologically 
linked to their initial resonant surface (see right snap- 
shot of fig. ©). At t = 478, the global reconnection 
process is completed and the islands disappear (step 

3) . Inspection of the maps of the electrostatic poten- 
tial just before the global reconnection (see fig. Q) 
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Fig. 8 Sx = tt/2. Structure of ^ ot (a;,t). (Top-Left) t 
380ta. (Top- Right) t = 390t a . (Bottom-Left) t 
430ta. (Bottom- Right) t = 440ta 
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Fig. 9 Sx = tt/2. Structure of Bl^(x,t). (Top-Left) t 
430ta. (Top- Right) t = 440ta. (Bottom-Left) t 
450ta- (Bottom- Right) t = 476r A 



shows that the presence of the m = mode, even if 
it appears to be much more complex, does not sub- 
stantially modify the plasma flow compared with the 
case when vq = 0, where typical butterfly-like struc- 
tures are easily identified. To clarify the role of the 
generated poloidal mean flow VQ at (x,t), it is interest- 
ing to focus on the time evolution of its structure. 
Fig. © shows VQ ot (x,t) at different times. The red 



lines represent vq — fQ Ot (a;,0) and the horizontal lines 
show the position of the resonances. Initially, we ob- 
serve that vq grows close to the resonances. Then, 
(step 1) occurs roughly when v < ' t (xi r ,t) has crossed 
zero (about t = 385) or some Alfven time after this 
has happened (about t — 395). Finally, the graphs 
t — 430Tyi and t = 4A0ta show that during the global 
reconnection, which starts at (step 2), the mean veloc- 
ity oscillates around zero in between the resonances, 
approaching, in some sense, the case without imposed 
shear flow. These observations show, first, that the 
nonlinear destabilization of the DTM occurs once the 
generated vq(x, t) compensates vq in the vicinity of 
the resonant surface, and second, that the global re- 
connection process starts and occurs once the mean 
velocity profile in between the islands oscillates around 
zero. 

It is instructive to follow the time evolution of 
the mean magnetic field B t ° t (x,t). Fig. © shows 
i?Q ot (x,t) at different times during the global recon- 
nection process, i.e. in between (step 2) and (step 
3). First, we can see that until (step 2), the pro- 
file and also the positions of the resonant surfaces 
(i?Q 0t (x,t) = 0) are roughly unchanged. This is fol- 
lowed by a flattening of the profile in between the 
resonant surfaces where the mean current decreases, 
while in the vicinity of the resistive layers, the pro- 
file becomes narrow signifying the generation of a 
strong mean current (see snapshot at t = AAOta and 
t = 450x^4, respectively). This effect is amplified un- 
til (step 3) occurs. It corresponds to the time when 
BQ ot (x,t) does not have a zero point at any location 
(see t = A76ta snapshot). In other words, there is no 
resonant surface anymore. 

Let us consider also the case Sx = 37r/4, keep- 
ing all the other parameters identical to those used in 
the case Sx — ir/2. The KH stable branch still ex- 
ists and shows the similar behaviour as the case with 
Sx = 7r/2. We find that the same behavior is observed. 
However, the situation is more complex because is- 
land rotation is initially involved in the dynamics. In 
fact the (step 1) event occurs at t = 660, when also 
the magnetic islands stop to rotate poloidally and the 
(step 2) event occurs at t = 1600, after a long pe- 
riod where the islands do not move poloidally (see 
fig. ( [TP)) ). Moreover, in this latter case, the global 
reconnection corresponds to a structure-driven non- 
linear instability of DTM characterized by an abrupt 
growth after a long-time-scale evolution [TJ] , such that, 
at t = 1610, the islands reach the wall. In both cases, 
however, the presence of the shear flow delays the time 
at which global reconnection occurs. The delay is of 
about 15% for Sx = tt/2 and 30% for Sx = 3tt/4. A 
much more systematic study should be performed, but 
this preliminary work shows that control of the global 
reconnection process through control of the shear flow 
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Fig. 10 Sx = 37r/4. Time evolution of the poloidal posi- 
tions of the center of the two islands 



in between the islands might be possible. 

6. Conclusion 

In this paper, we have presented a study of the 
role of a shear flow in between magnetic islands on 
the global magnetic reconnection processes. We first 
have identified the different linear regimes which exist 
when a shear flow is present in between two surfaces 
where DTM develop. We have seen that according to 
the distance between the resonant surfaces and the re- 
sistivity of the plasma different regimes exist. When 
the distance is larger than the typical magnetic shear 
length, no global reconnection occurs and the islands 
rotate 15 . In this paper, we did not focus on those 
regimes. When it is of the order or smaller than the 
typical shear length of the system, an island can be 
locked or rotate linearly according to the value of Sx. 
Moreover, the system can be KH stable or unstable 
according to the resistivity and, of course, to the am- 
plitude of the shear flow. We have shown that, inde- 
pendent of the regime (rotating island regime or not) , 
a breaking of symmetry linked to the presence of the 
shear flow leads to the generation of a mean poloidal 
flow at early times. We have shown that flow devel- 
ops first in the vicinity of the resonant surfaces and 
is at the origin of the radial displacement and desta- 
bilization of the DTM structure. We have found that 
the presence of the shear flow delays the global re- 
connection process, suggesting that its control might 
be usefull for controlling the DTM. Further investiga- 
tions are, however, necessary to precisely evaluate the 
impact of the shear flow and its amplitude, and also 
the impact of the KH instability in the low resistivity 
regimes, on the global reconnection processes of the 
double tearing mode. 
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